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Abstract 

We present a model for high-energy solar flares to explain prompt proton and electron acceleration, which occurs 
around moving X-point magnetic field during the implosion phase of the current sheet. We derive the 
electromagnetic fields during the strong implosion of the current sheet, which is driven by the converging flow 
toward the center of the magnetic arcade. We investigated test particle motion in the strong electromagnetic fields 
derived from the MHD equations. It is shown that both protons and electrons can be promptly (within 1 s) 
accelerated to — 70 MeV and ~ 200 MeV, respectively. This acceleration mechanism can be applicable for the 
impulsive phase of the gradual gamma ray and proton flares(gradual GR/P flare), which have been called two- 
ribbon flares. 


1. Introduction 

One of the important problems in high-energy solar flare physics is to understand the physical mechanism(s) of high 
energy particle acceleration during solar flares. From the observations by SMM and Hinotori satellites, it became 
clear that there are two classes of gamma-ray/proton(GR/P) flares- impulsive GRIP flares and gradual GRIP flares( 
see, for a review , Bai and Sturrock 1989). The GRIP flares refer to flares that produce nuclear gamma-rays and /or 
energetic interplanetary protons. Most short flares(< 100 s) corresponds to be impulsive GRIP flares, while most of 
the long-duration flares(> 200 s) corresponds to the gradual GRIP flares. In the impulsive GR/P flares there are two 
phases of particle acceleration : the first phase is that both electrons and protons are accelerated to — 10 MeV and to 
- 100 MeV energies , respectively, within 1 s (Kane et al. 1986; Riger 1989), and second phase is that electrons up 
to 100 MeV and protons up to GeV energies are accelerated within a few seconds to 100 s. 

In order to explain the above prompt particle acceleration to relativistic energies, Sakai(1990) showed that during 
the 3D X-type current loop coalescence, where two crossed flux tubes interact in one point, both protons and 
electrons can be promptly( within less than 1 s) accelerated to - 100 GeV and ~ 100 MeV, respectively. De Jager and 
Sakai (1991) showed that the burst duration of elementary flare bursts(5-25 s) observed during the impulsive GR/P 
flares can be explained quantitativ ely by the mechanism of the 3D X-type current loop coalescence. 

In the present paper we show a model for the long-duration gamma ray /proton flare(the gradual GR/P flares) to 
explain prompt proton and electron acceleration during the impulsive phase. These flares have been called two- 
ribbon flares and almost all gradual GR/P flares show impulsive behaviour in the beginning of the hard X-ray 
emission and gradual behaviour later on (Bai and Sturrock 1989 ). In most large two-ribbon flares with filament 
eruptions the filament begins to move up several minutes or more before the onset of the impulsive phase(Martin 
and Ramsey 1972; Kahler et al. 1988). Recently Sakai and Koide (1991) presented a theory of the filament eruption 
before the impulsive phase of solar flares. Figure 1 shows a schematic picture of the global magnetic field 
configuration and the filament which is located above the X-point magnetic configuration near the center of the 
magnetic arcade. Sakai and Koide (1991) showed that the upward motion of the X-point which may trace the 
filament eruption begins several minutes before the impulsive phase where the magnetic field, Bx produced by the 
current near the X-point can exceed the magnetic field, Bz along the filament. If the Bx exceeds the Bz, then the 
explosive magnetic reconnection (Tajima 1982; Sakai et al. 1984; Sakai and Ohsawa 1987) can be triggered in the 
current sheet. During the implosion phase of the current sheet, the strong electric field can be induced by the strong 
converging flow around the X-point. 

In Section 2 we drive the electromagnetic fields during the implosion phase of the current sheet from the MFTD 
equation. In Section 3 we investigate a test proton and electron motion under the electromagnetic fields derived in 
the previous section. We will show that both protons and electrons can be promptly accelerated within 1 s to ~ 70 
MeV and - 200 MeV, respectively. In Section 4 we summarize our results. 
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2. Electromagnetic Fields around Moving X-point 


2.1 MHD Equation and Energy Equation 

In this section we present basic equations describing the dynamics of the current sheet near the magnetic X-point, 
which is located under the filament as shown in Figure 1. We assume that the filament is supported in die upper part 
of the X-type magnetic field along the z-direction, in which there exists the magnetic field, Bz. Following the 
gradual shearing motion at the footpoint of the magnetic arcade, the plasma upward motion can be induced, and e 
current near the X-point can be increased due to the converging flow near the current sheet The current sheet 
continues to elongate to the vertical direction, and to become thin in the horizontaJ direction. We will focus on die 
dynamics of the plasma around the X-point to determine the electromagnetic fields during thinmng of ‘hecurrent 
sheet. The x-axis is taken in the vertical direction, and the y-axis is in the horizontal direction. The current sheet 
assumed to be homogeneous in the horizontal z-direction. ... rr 

We begin from the following MHD equations including the gravity and Joule heating effect in the energy equation. 


| + v .(pv)-° 


pfil + v-Vv! = -VP+ -— ( V x B) x B - pge x 

dt ) 

— =Vx(vxB)+-^-—AB 

dt 4 no 

dP (V x B) 2 

— + (vV)/> + yPV-v-(y-l) i — — 

dt 16 no. 


( 2 ) 

( 3 ) 

( 4 ) 


where r , v, P, B are the density, velocity, pressure, and magnetic field, respectively, and y is * e adiabatic constant. 
The gravitational acceleration is given by 
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where R = R s + r 0 and go = GM S R 2 . M s , R s and r 0 are the solar mass, solar radius and distance to the 
origin of the coordinate from the solar surface, respectively. The right-hand side in Equation (4) shows the effect of 
the Joule heating. 


2.2 Electromagnetic Fields around the Current Sheet 

During the gradual shearing motion of the photosphere, the current near the magnetic X-point can increase due to 
horizontal plasma inflows from both sides. This horizontal plasma flow, v y around the X-pomt can be approximately 

given by 


where a(t) is a time dependent scale factor. The scale factor a(t), which is determined later, characterizes continuous 
change of thickness of the current sheet due to the plasma inflows. The vertical flow, v x is taken to be 
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where vxO(t) and another scale factor b(t) are determined self-consistently later. 
The components of the magnetic field are taken as 
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B z = B l0 (t) 


( 10 ) 


where k is a characteristic scale-length of the current sheet, which can be considered to be the same as the radius 
of the filament located above the X-point. The magnetic field B n in Equation (9) shows a role that the X-point can 
move during the evolution of the current sheet. Unknown functions, B x o, B n> ByO and B z o can be determined, self- 
consistently later. 

Substituting Equations (6) and (7) into Equation (1), we find that the density p(t) is only a function of time and 
given by 


P (t) = — &— 
<<Mt) 


(ii) 


where po is a constant. 

From the induction equation (3), using the expression for the magnetic fields (8)-(10) and velocities (6) and (7) > 
we find 
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where B 0 and Bqq are constants. 

We assume the pressure P(x,y,t) to be 
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We here summarize normalized basic equations , which will be numerically solved in the next section. 
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where Pp = c s 2l.ap * is plasma be.a-ra.io, which is determined from the magnetic field. Bo * £ 

current induced near the X-point. The magneuc Reynolds number, S is S- Tr/Ta ( XB ) 

parameters Gi in Equation (18) and G 2 in Equation (20) are given as 
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In 


the above equations (17)-(22), the time , B n 0, v x 0. % POO are normalized by x A = X / vap , Bo. va P , Pa and Pa, 


The normalized time-dependent part of the magnetic fields, Bx and Bz are given by 


1 

BjU “ 2 

a 


(25) 


B « = R >V b 


(26) 


where Rb= Boo/Bo gives the initial ratio between Bz and By ... . 

m deTc field E induced by the change of the magnetic field B( eqs. [8]-[10] ) can be determined from 


E - -v x B / c 


(27) 


By use of equations (6)-(10), we obtain the components of the inductive electric field; 
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We can obtain numerical results of the electromagnetic fields by solving the basic equations (17)-(20). The initial 
conditions are taken as a = 1, da dt = -0.001, b = 2 , db/dt =0, Bu = 0.01, v x q = 0.001, and other parameters are 

taken as Gi = 0.01, G 2 = 0.27, S = 10 5 , and Rb = Boo/Bo = 100/10. The plasma beta-ratio p p is 1.0 and the 
adiabatic ratio is y = 5/3. 

Some components of the inductive electric fields become explosively strong near the time 1.18 ta The reason is 
due to the fact that near the implosion phase, the converging flow, v y is increased by about 100 times and the 
magnetic field, B xa is also increased by about 150 times through the current pinch. Therefore, the inductive electric 

field E x = v v B x () c can be amplified to 1.5 x 10^ times, compared with the initial values. This inductive electric 
field before the implosion is about 70 V/m , if we take v ap = 7x 10 6 cm/s and B 0 = 10 G. The inductive electric 
field near the implosion can become E x ~ 10 6 V/m. This strong electric field near the implosion is important for 
the high-energy particle acceleration. If we take X. = 2x1 0 9 cm, the Alfven transit time t a is t a = A/v ap = 300 s. 
Then the duration of implosion is about 0.01 t a = 3 s. As seen in the next section, the acceleration time for both 
protons and electrons is quite rapid compared with the duration of implosion , 3 s. Therefore we will investigate the 
motion of a test particle under the electromagnetic fields (28)-(30) in the next section. 


3. Proton and Electron Acceleration 


3.1 Equation of Motion 

We consider the motion of a test proton under the electromagnetic fields given by equations (28)-(30) 
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whereP =P/m p c, E =E/Eo, B = B/ Bo, R p = c/ctfcpk, x = x / k, y = y / X , z= z/ k. 

The time is normalized by the proton cyclotron period co _1 cp (T- to cp t , to cp = efio m p c ) 
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3.2 Numerical Results for Proton Acceleration 

The linie-dependenl amplitude* of the electromagnetic field* are assumed to be instant during 

r5r^7.^:^ 

,3 r hz, o.o, 9) - - L.n :: 

initial conditions with x(0)=y(0)_z(0)_U.5, rx-ty- 1 z ll £ . T^nmioncanbe 

proportional to the time. Therefore we can evaluate the time dependency of the momentum P x as 
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for y(0) = 0.5. Similarly, we find 
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for y(0) = 0.1 and 


P x = 4.0 x 10^(0^/ + 10" J 
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implosion time the proton can be accelerated to ~ 70 MeV for protons with >(0)> 0.05. 

3.3 Numerical Results for Electron Acceleration 

In this section we present numerical results for electron acceleration by solving the following equations of motion 

for a test electron. , ^ , * , v . R _ io-8 

The time is normalized by the electron cyclotron period co" ce ( t - Wee t . % ! - ' «®0 ** > J^-4 The 

and v Ic - H 4300 The initial conditions for the electron are x(0)-y(0)-z(0) 0.., x l y z. 

dependency of the average momentum P x as 
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for y(0) = 0.5, 


p = 1.0 X 10' S «U + 10" 


for y(0) = 0. 1 and 


= 4.0 x 10 + 10 
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ce , V ’ A 4 .1 

for y(0) = 0.05. These results are almost same for the initial momentum in the range of P(0) 1 “ ^ 
as for the other initial positrons, x(0) and z.(0). As the electron cyclotron frequency ts o>ce = H* s for B 0 - 1 . 
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we conclude from Equation (38) that within 1 s the electron with y(0)> 0.05 can be accelerated up to ~ 200 MeV. 


4. Summary 

In this paper, we presented a model for high-energy solar flares to explain prompt proton and electron acceleration 
which occurs around moving X-point magnetic field. We derived the non-stationary electromagnetic fields during 
the strong implosion of the current sheet, which is driven by the converging flow toward the center of the magnetic 
arcade .We investigated a test particle motion in the electromagnetic fields derived from the MHD equations. It has 
been shown that both protons and electrons can be promptly (within 1 s ) accelerated up to ~ 70 MeV and ~ 200 
MeV, respectively. This acceleration process can be applicable for the impulsive phase of the gradual gamma ray 
and proton fiares(gradual GR/P flare), which have been called two-ribbon flares. 
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Fig.2 Time history of the proton's Lorentz factor. 


Fig. 1 Schematic magnetic configuration and coordinate system. 
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Fig.4 Time history of the electron's Lorentz 
factor( upper) and momentum P x (lower). 
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Fig.3 Time history of the proton's 
momentum. 
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